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The localization lengths A of one-dimensional disordered systems are studied for electronic wave- 
functions in the Anderson model and for vibrational states. In the first case, the site energies e and in 
the second case, the fluctuations of the vibrating masses m at distance £ from each other are long-range 
correlated and described by a correlation function C(£) ~ £~~" with < 7 < 1. In the Anderson 
model, we focus on a scaling theory that applies close to the band edge, i.e. at energies E close to 
2. We show that A can be written as A = Xof-,(x), with A = ( e 2 }- 1 /(4-7) ^ X (E = 2, (e 2 )), 



= A (2 — E) and the scaling function fj(x) = const for 1 < 1 and / 7 (a 



,(3-7)/2 



for 



i>1. Mapping the Anderson model onto the vibrational problem, we derive the vibrational local- 



ization lengths for small eigenfrequencies u), A ~ (m) 



(3-7)/2 



"(1+7) 



where (m) is the 



mean mass and (m 2 ) the variance of the masses. This implies that, unexpectateley, at small u, A is 
larger for uncorrected than for correlated chains. 



L INTRODUCTION AND MODELS 

A large amount of work has been done in the past 
decades to understand localization behavior in randomly 
disordered chains. Most work has been concentrated 
on uncorrelated [[[]] and short-range correlated disorder 
( see [||] and references therein). In this work, the effect 
of power-law long-range correlated disorder on the lo- 
calization properties of linear electronic and vibrational 
systems is discussed. 

First we consider the Schrodinger equation of a 
single-particle electronic wavefunction in tight binding 
approximation (Anderson model) [^J, 



1p n+ l + 1p n - 



-(2-E + e n )i) n . 



(1) 



Here E is the energy and \tp n | 2 is the probability to find 
an electron at site n, with n € {1, AT}. The e„ are the 
site potentials and are random numbers, taken from an 
interval of width A, e„ € [-A/2, A/2] (diagonal dis- 
order). (e 2 ) = (1/N) Y^,n=i e n i s their variance and the 
average value (e) = (1/N) J2n=i e ™ ^ s zero - Th e term 
— 2ip n on both sides is introduced for technical reasons. 

Second, we investigate the related problem of the 
localization length of vibrational states in disordered 
chains. In this case, the wavefunction ip n of Eq. ([!]) is 
substituted by the vibrational amplitudes u n of N par- 
ticles at sites n, coupled by unit spring constants be- 
tween them. This problem is described by the time- 
independent eigenvalue equation with the vibrational 
amplitudes u n and the eigenfrequencies u>: 



-1 - 2u n = -m n uj u n 



(2) 



Contrary to the e n of Eq. ([[]), all masses m„ are positive 
and taken from an interval m n 6 [(m) — A/2, (m) + 
A/2], where (m) = (1/N) Yln=i m » > is their mean 
(l/^)E»=i(»»n - M) 2 the vari- 



value and (m ) 
ance. 



The localization length A in the Anderson case is 
defined by the exponential decay of the wavefunctions 
4> n , lim„^oo I V'n/V'o I = exp[— n/A]. In the vibrational 
problem, tfju/fpo is substituted by u n /uo, which decays 
accordingly. 

In the standard Anderson model with uncorrelated di- 
agonal disorder it has been recognized since long that in 
d = 1 all states are exponentially localized i.e., 
that A approaches a constant for large N, A weak disor- 
der expansion yields A = (e 2 }- 1 / 3 /((2 - E)/(e 2 ) 2 ^) 
close to the band edge while at the band cen- 

ter, a Green's function technique yields A ~ (e 2 ) 1 
[pp. In some distance from the band center, the diag- 
onal elements of the Green's functions can be devel- 
oped by a second-orderperturbation theory, yielding 
X(E) ~ (4 - £ 2 )/<e 2 } [§. 

In the related vibrational problem with uncorrelated 
random masses, it is well accepted that all states in 
d = 1 are localized, except for the case uj = 0, where 
the disorder term m„w 2 disappears. Transfer matrix cal- 



culations in d = 1 yield A 



for small lo [1( 



Here, we concentrate on long-range correlated disor- 
der, characterized by a correlation function C(£) ~ £~~< . 
In the Anderson case, C(£) describes the correlations 
between site energies e n and e n +i at sites n and n + t, 



1 N 

C(l) = lim — e n e n+e = (e n e n+ i) 

N—>oc 1\ — ' 



(3) 



with the correlation exponent 7, < 7 < 1. In the vi- 
brational problem, the mass fluctuations are long-range 
correlated and described by C(l) = (m n rh n+ i) = £ 1 
with m n = m n — (m). The case 7 > 1 describes 
only short-range correlated potentials. From random 
walk theory, we expect that series with 7 > 1 fall into 
the same universality class as uncorrelated potentials 



For the one-dimensional Anderson model of Eq. ([]]) 
with long-range correlated potentials, a scaling form for 
the localization length X(E, (e 2 )) close to the band edge 
E = 2 has been developed recently and supported by 
preliminary numerical simulations [[13]]. For the vibra- 
tional problem, the localization length in the presence 
of correlations has not been adressed yet. In this paper, 
we determine explicitely the scaling function for the An- 
derson system and derive its asymptotic form by scaling 
arguments (section II). This leads also to an explicit ex- 
pression for A at intermediate energies in some distance 
from the band edge (cf. Eq.(|lC|)). In section III, we map 
the Anderson model onto the vibrational problem and 
show that for small u>, A(w) ~ w~' 1+7 ^ which agrees 
with the result A ~ lu~ 2 of [|l0|] for random uncorre- 
cted potentials, if we describe them by 7 = 1. Accord- 
ingly, contrary to the expectations, A is decreased by the 
presence of long-range correlations for u> — > 0. 

II. THE LONG-RANGE CORRELATED ANDERSON 
MODEL 

In a previous work [^3j], it was shown by a space- 
renormalization procedure that at the band edge, E = 2, 
the localization length A scales as 

A(£ = 2, (e 2 )) ~ (e 2 )" 1 /^) ^ Ao; (4) 

For E < 2, but still in the neighborhood of 2, A depends 
on both, Ao and 2 — E and can be written as: 

X(E,(e 2 )) = A / 7 (z), (5) 
where x = A 2 , (2 - E) = (2 - E){e 2 )- 2 / {i -^ (6) 

and f 1 (x) is a correction function that depends on 7 and 
approaches a constant for x — ► 0. For 7=1, one recov- 
ers the exponent 1/(4 — 7) = l/3> in accordance with 
the result of [0] for uncorrected potentials. 

We are interested in the behavior of / 7 (x) in the limits 
of large and small values of x. For uncorrected poten- 
tials e n , we know already that / U ncorr(^) — * const for 
x < 1 @ and / unC orr(aO -^xforx»l [§|]. We 
therefore assume a power-law behavior 

/ 7 (a;) ~ x a (7) 

for correlated potentials, where a depends on 7. 

To derive the exponent a, it is helpful to realize that 
1 / (2—E) is proportional to the square of the wavelength 
A of the electronic wavefunction ((e) = 0), 

A 2 ~l/(2-£). (8) 

This can be seen most easily in an ordered chain, where 
all e n = and Eq. (|J) reduces to ijj" (x)/ip(x) = 
— (2 — E). The solution is a harmonic function with the 
wavelength A. In the case of disorder, we have calcu- 
lated many states of correlated and uncorrected chains 



in different disorder intervals A by the Lanczos algo- 
rithm and have verified by Fourier transformation that 
relation ([$]) still holds. So, A describes the oscillating 
part of the wavefunction, while A describes the expo- 
nential decay of the envelope. A similar relation for the 
wavelength has been found for vibrations in percolation 
clusters in d = 2 and d = 3 in [|l4]] . 

Next, we look at a wavefunction with a given wave- 
length A and discuss, how A is influenced by a change 
of the variance of the potentials (e 2 ), i.e. by a change of 
Ao- We discuss the cases i«l and x 3> 1 separately, 
which correspond to A ^ Ao and A <C Ao, respectively. 
We have to keep in mind that A is solely determined by 
E, whereas Ao is determined by (e 2 ) and 7. 

(i) For A 3> Ao, close to the band edge, the wave am- 
plitude decays within the first wavelength A. Accord- 
ingly, the wavefunction is independent of A and would 
not change for A — > 00, i.e., by approaching the band 
edge, where A ~ Ao = (e 2 ) -1 /' 4 " 7 ). Hence we ex- 
pect for all 7 that the scaling function / 7 (x) = X/Xq 
approaches a constant for small x. 

(ii) In the opposite case, A <C Ao, the wave amplitude 
performs many oscillations before the envelope function 
decays completely. It is reasonable to assume that for 
the dependence of A on the potential landscape, long- 
range correlations are not relevant, when the wavefunc- 
tion oscillates rapidly. In this case, the term e n ip n on 
the right-hand side of Eq. ([!]) oscillates rapidly for both, 
correlated and uncorrected potentials. So, in the case 
A <C Ao, we expect the same dependence of A on (e 2 ) 
for correlated and for uncorrected potentials. For the 
latter, we know that in some distance from the band 
edge, the behavior of A must cross over to the result 
°f [f^L A ~ (e 2 ) -1 and in view of the preceeding re- 
mark, we expect the same behavior also for correlated 
potentials. Inserting this relation into (|J) and (Q), we 
find f 1 (x) ~ x a , with a = (3 — r y)/2. In summary, 
we find for the asymptotic cases of the scaling function 
fry (x) the power-law behavior 

:>>: <» 

Inserting Eq. into (|J) and (^), we find for the en- 
ergy dependence in the case A <C Ao with the additional 
constraint 2-E<l 

X(E, (e 2 )) - (e 2 )- 1 (2 - E)^)/ 2 . (10) 

Accordingly, we can distinguish between three energy 
regimes, (a) For 2 — E <C (e 2 ) _2 /( 4 >\ i.e., very close 
to the band edge, we have A ~ (e 2 ) _1 /( 4-7 ), indepen- 
dent of E. (b)Forl > 2-E > (e 2 )- 2 /^\ Eq. @ 
applies, (c) For still larger values of 2 — E, the scaling 
behavior breaks down. 

For growing values of 2 — E, it was obtained that in 
some distance from the band center A (E) ~ (2 — E)(2 + 
E)j (e 2 ) for uncorrected potentials |^]. This expression 
emerges into Eq. (|To|), if we set 7 = 1, which yields 
A ~ (2 — E)/(e 2 ). So we can see again that the case of 
uncorrected potentials is well described by 7 = 1. 



E -> w 2 (m) 



and 



w 2 m„. (11) 




o (2-E) 



FIG. 1.: Plot of the scaling function / 7 (a;) = \(E)/\o versus 
a; = A 2 , (2 — _B) close to band edge for the Anderson model 
in d — 1. The curves show (from bottom to top): uncorre- 
cted disorder (large open symbols), correlated site potentials 
of 7 = 0.5 (large filled symbols), 7 = 0.1 (small open sym- 
bols) and 7 = 0.01 (small filled symbols). The symbols stand 
for 2 - E = 10~ 4 (circles), 10~ 3 (squares), 10" 2 (diamonds), 
0.1 (triangles up), 0.2 (triangles down), 0.3 (x) and (e 2 ) was 
varied between 10~ 5 and 10 _1 . The typical average of A was 
performed over 1000 chains of length N = 2 20 . The straight 
lines of the theoretical slopes (3 — 7)/2 are guides to the eye. 

In order to test this scaling function, we have gener- 
ated v — 1000 linear chains of size TV = 2 20 by the 
method of Fourier transformation [ jl5| ] and calculated 
the individual localization lengths \ by the transfer- 
matrix method ji|, p^|[r^ ] . To obtain the mean localiza- 
tion length, we performed the typical (logarithmic) av- 
erage X(E, (e 2 )) = A typ = exp[i Yh=i ln A *]- 

To obtain the scaling function / 7 (x), we plotted A/Ao 
versus x = A 2 , (2 — E). The results are shown in 
Fig. 1 for (from bottom to top) uncorrelated potentials 
(described by 7 = 1) and correlated potentials with 
7 = 0.5, 0.1 and 0.01. Five values of (2 — E) between 
10~ 4 and 0.3 were calculated while (e 2 ) varied between 
(e 2 ) = 10~ 5 and 10 _1 . The figure supports the scaling 
theory. The data for the different 2 — E and (e 2 ) fall 
onto the same scaling functions f 1 (x), which show the 
expected behavior (cf. Eq. (^|)) in the asymptotic cases. 
The lines are guides to the eye with the theoretical slopes 
of a = (3 — 7)/2. A crossover close to x — 1 separates 
the cases x < 1 and x > 1, i.e. A 2 , < 1/(2 - E) 



and A 2 , > 1/(2 - E). In the case x < 1, we find 
A/Ao = const, i.e. a = 0. 



III. THE LONG-RANGE CORRELATED 
VIBRATIONAL MODEL 



This scaling theory can be generalized to the vibra- 
tional problem, described by Eq. (Q). We split m„ 
into the average part (m) and the fluctuation part fh n , 
m-n — (m) + rh n and can map the Anderson problem 
onto the vibrational problem by replacing 



We obtain the wavelength A ~ ((m) 1 / 2 uj)~ x by insert- 
ing ( p"lj ) into (Q) and the scaling variable x by inserting 
([□]) into (|[), yielding (with (m 2 ) = (m 2 )) 

x= <™> W - 2 ^ 4 -^<m 2 r 2/(4_7) - (12) 
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FIG. 2.: The scaled localization length A (m 2 ) (m)" (3_7)/2 
for the vibrational chain is plotted versus the eigenfrequency 
uj for uncorrelated masses (circles) and correlated masses of 
7 = 0.1 (diamonds) for different (m) and disorder intervals 
A: (to) = 0.5, A = 0.99 (large open symbols), (to) = 0.5, 
A = 0.5 (filled symbols) and (to) = 0.25, A = 0.49 (small 
open symbols). The values of A = 0.99 and 0.5 refer zu 
(m 2 ) ~ 0.008 and 0.002, respectiveley. The typical average 
of A was performed over 1000 chains of length N = 2 23 . The 
lines with the slopes 1 + 7 are guides to the eye. Small devia- 
tions from the straight line at very large values of A ~ N/2 are 
finite-size effects. 

Since (m 2 ) < (m) 2 , we can easily verify that for 
u 2 < l/(m), the scaling variable x is greater than 1. 
Therefore, contrary to the Anderson problem, we have 
only one scaling regime x > 1. Inserting ([FT]) into (|lO[), 
we find A for the vibrational problem: 



(13) 



This relation is valid for small values of Lo 2 (m). In our 
simulations, we found that for e.g. (m) = 0.5, it holds 
approximately for to < 1. In order to test Eq. (JT^), we 
have calculated the typical average At yp of 1000 chains 
of length 2 23 for different values of 7, (m) and (m?) 
by the transfer-matrix method. In Fig. 2, the scaled lo- 
calization length A (m 2 ) (m)~( 3 ~ 7 )/ 2 is plotted versus 
oj for 7 = 0.1 and for uncorrelated random chains (de- 
scribed by 7 = 1). Both, the data collapse for chains of 
different (m) and (m 2 ) and the slopes support Eq. (pj|). 
The lines are guides to the eye with the theoretical slopes 
of —(1 +7) and describe the behavior of A very well. 

Two phenomena are worth mentioning here. First, A 
diverges if lo — > 0. The reason is that the disorder term 
m n in (ph appears only in combination with the eigen- 
value uj and therefore disappears for ui — > 0. Second, 



for large ui, A is larger for the correlated than for the 
uncorrelated chains. However, below some crossover 
frequency (that depends on 7), this behavior changes 
and the vibrational excitations of the correlated chains 
become more localized. This, at the first glance unex- 
pected behavior can be understood as follows. 

Large oj correspond to small wavelengths A. Particles 
m n in distances less than A/2 from each other move 
into the same direction and can roughly be considered 
as an effective hyperparticle Mj ~ X^n=i TO ™' where 
the sum runs over all /i particles in the region of the size 
of A/2. The effective disorder seen by the vibrational 
excitation depends on the distribution of the Mj in the 
region, where the wave amplitude is large. For strongly 
correlated chains, there exist large regions of masses m n 
that are either below the average mass (m) or above. So, 
as long as A/2 is smaller than the size of these regions, 
the distribution of the Mj is more narrow for correlated 
than for uncorrelated chains, which leads to larger A. 

This behavior changes when the size of A/2 is well 
above the size of the correlated regions. Now, on scales 
of A/2 the system looks much more disordered than be- 
fore. We know from random-walk theory that for large 
H the variance of the Mj scales as ~ /i 2 ~ 7 Hence, 
for large /1 the distribution Mj is broader in correlated 
than in uncorrelated chains and therefore, the disorder 
seen by the vibrational excitation is even larger. There- 
fore, we expect smaller A for the correlated chains, when 
u -> 0. 

The same effect appears also in the Anderson model, 
where, according to Eq. (10), A is also smaller for cor- 
related than for uncorrelated chains in the limit of small 
values of 2 — E. For systems at the band edge, this has 
already been reported in [O]. 



IV. CONCLUSIONS 

In summary, we have studied the effect of long-range 
correlated disorder on the localization lengths of elec- 
tron states in the one-dimensional Anderson model and 
of the vibrational states of harmonic chains. In the 
Anderson model, the site potentials e are long-range 
correlated, while in the vibrational problem we have 
long-range correlations in the fluctuation part m„ of the 
masses around their mean value (m). 

First, we studied the scaling behavior of the localiza- 
tion length A in the Anderson model close to the band 
edge E — 2. We have shown that two characteristic 
lengths occur: the wavelength A that describes the os- 
cillating part of the wavefunction and the localization 
length Ao that describes how the envelope of the wave- 
function decays at the band edge. While A depends on 
the energy E, Ao depends on the variance (e 2 ) and the 
correlation exponent 7. We developed a scaling theory, 
with the scaling variable x — A 2 , /A 2 . Using qualitative 
arguments, we derived the asymptotic behavior of the 



scaling function, both, for A ^ Ao and A < Ao, and 
thus the behavior of A in these limits. It was found that 
A is independent of E at very small values of 2 — E. For 
larger distances from the band edge, but E still in the 
neighborhood of 2, it depends on 2 — E and on (e 2 ) by 
power-laws. 

Second, we developed an analogous scaling theory 
for the vibrational problem by mapping the Anderson 
equation onto the equations of motion of harmonic vi- 
brations. We found that X(uj) ~ tj - ( 1+7 ) in the limit of 
small frequencies and for < 7 < 1. This relation may 
constitute a possibility to measure correlations by deter- 
mining the localization length of vibrations. It implies 
that A grows faster for uncorrelated than for correlated 
chains, when u> — > 0. This is in close analogy to the be- 
havior of A in the Anderson model at the band edge and 
arises from the coupling of correlated masses over very 
large distances, which leads to larger fluctuations than 
in the case of uncorrelated potentials. Accordingly, for 
very small frequencies, correlated chains have smaller 
localization lengths than uncorrelated ones. 
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